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Abstract. We consider an interacting particle system on a graph which, from a macroscopic point 
of view, looks like 7L d and, at a microscopic level, is a complete graph of degree N (called a patch). 
There are two birth rates: an inter-patch one A and an intra-patch one <j>. Once a site is occupied, 
there is no breeding from outside the patch and the probability c(i) of success of an intra-patch 
breeding decreases with the size i of the population in the site. We prove the existence of a critical 
value \ C r(<t>, c, N) and a critical value C r(A, c, N). We consider a sequence of processes generated by 
the families of control functions {ci}i 6 N and degrees {iVijigN; we prove, under mild assumptions, 
the existence of a critical value i cr (\,<t>,c). Roughly speaking we show that, in the limit, these 
f^J • processes behave as the branching random walk on l d with intra-neighbor birth rate A and on-site 

birth rate <j>. Some examples of models that can be seen as particular cases are given. 

Keywords: contact process, restrained branching random walk, epidemic model, phase transition, 
critical parameters. 

AMS subject classification: 60K35. 
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q\ • 1. Introduction 

Stochastic models for the demographic expansion (or contraction) of a biological population have 
attracted the attention of many researchers since Galton and Watson introduced the branching 
process as a very simple and non-spatial model to study the survival of surnames. After the 

oo ; 

appearance of interacting particle systems and of the paper of Harris [10] on the contact process 
j> , it has become increasingly clear that the spatial structure of the environment and the interaction 

•rH . 

^ between individuals bring a desirable complexity in the models one can define (see for instance [9] 
■ for a discussion on spatial models). 

We want to introduce and study a model for a biological population where individuals breed and 
die living on a patchy habitat. Clearly, if we substitute the words infect and recover for breed and 
die, our model also serves as a picture for the spreading of an infectious disease. 

The main idea is that the environment is subdivided in patches and the patches are centered 
at the vertices of the d-dimensional lattice Z d (see Figure Q3 where each square is a patch and 
d = 2). Individuals in the same patch are strongly connected to each other (indeed the patch 
may be thought as a complete graph with an intra-patch breeding rate 4>) and are also related to 
individuals in neighboring patches (the inter-patch breeding rate being A). 

Each patch has the same fixed number of sites and a site accepts individuals from neighboring 
patches if it is empty. This kind of biparametric interaction is similar to the one introduced in 
[TT] and later modified in 0. In both papers the population is divided in equally sized clusters, 
one cluster per vertex of Z rf ; in [11] individuals in the same clusters die at the same time, while 




a patch of 
N=9 sites 



Figure 1. The habitat. 



in [3] they die one at the time. In their model the cluster has no inner structure, that is, the 
only important thing is whether it is empty or not (if not, immigration is not allowed). In our 
model instead of clusters, we have patches and the process depends not only on the total number 
of particles in the patch, but also on their position among the sites within the patch. For instance 
let there be n particles in the patch. If they are on the same site, then all sites in the patch 
but this accept immigration (that is, reproductions from neighboring patches) while intra-patch 
reproductions are, a priori, always allowed. If they are scattered in n different sites all this sites 
will block immigrants. 

To strengthen the dependence of the process on the number of individuals per site, we admit a 
self-regulating mechanism, namely there exists a control function c £ [0, 1] N such that if a site is 
already occupied by j individuals, newborns are accepted only with probability c(j). This kind of 
regulating mechanism has been introduced in [5] as a general model of which the branching random 
walk and the contact process are particular cases (when c = 1 and c = Sq respectively) . Note that 
the presence of this function c mimics the fact that when an ecosystem is already exploited by 
many individuals, the younglings may not have enough food to live and to reach sexual maturity, 
thus many of them will not contribute to the breeding process. On the other hand many animals 
feed on their younglings if food is scarce. 

We observe that the idea of a patchy environment has already been introduced in [3] . The main 
difference is that in [JJ the inter-patch interaction is allowed only between the centers of neighboring 
patches, while in our model each site may receive immigrants from the neighboring patches. 

In the present work we mainly study the effect of the parameters on the probability of survival 
of the entire population. In particular we proved if A and <p are both small then there is almost 
sure extinction (Theorem I3.ip . This result is in agreement with the conclusion of the mean field 
model in a particular case (see Subsection I2.4p . If is small then survival is guaranteed only by A 
sufficiently large; while <j> large and a "mild" control function c implies survival even when A = 
(the precise statements are Theorem 13.21 and Corollary 13.31 where we give explicit conditions on (f> 
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and c). Analogously a large A guarantees survival even when (f) = 0, while for small A a large <f> 
(and a suitable c) is needed to ensure survival (see Theorem I3.4|) . 

To investigate the impact of different environment structures on the process we let the patch 
dimension go to infinity. We obtain Theorem I3.5f l) which agrees with the intuition that the 
limiting process is a branching random walk with intra-site breeding rate eft and inter-site rate A. 
The reader may want to compare this convergence to the one studied in [8], where a long range 
contact process is proved to converge to a branching random walk. We also study the effect of a 
progressive removing of the control function: we obtain survival when <ft > 1 and c is close to 1 
(Theorem 13.51 (2) ) . Finally we study the asymptotic behavior of the critical parameter of the general 
process and of some particular cases when the size of the patch goes to infinity (Proposition 13.71 
and Corollary 13. 9p or when the control function is progressively removed (Proposition 13. 8| ). 

Here is a brief outline of the paper. In Section [2] we introduce our process and discuss some 
models that can be seen as particular cases. Subsection 12.41 is devoted to the study of the mean 
field model. In Section[3]we discuss our main results; the proofs can be found in Section^ Section[5] 
is devoted to the discussion of open questions. 

2. The model 

The patches of the habitat are complete graphs K/y with iV vertices (that is any vertex is a 
neighbor of any other vertex, including itself) . In each single patch the sites are enumerated from 1 
to N and in each site there might be a priori any number of particles. The patches are located at the 
vertices of the lattice Z d . Formally, the state space is N £ n where the environment is £fr := Z rf x K/v 
and the graph structure of £fj is as follows: (x, r) ~ (y, r\) if and only if \x — y\ Z d < 1 (where | • | Z d 
denotes the usual distance in the graph Z rf ). On the other hand, if x,y £ Z rf then we write x ~ y 
meaning that \x — y\ Z d = 1. With this notation, the patch at x £ Z d is the set {(x,r) : r £ K/v}. 

The dynamics of the process is the following. Each particle dies at rate 1; when a site is vacant 
it receives births from each particle in the same patch at rate 4>/N and from each particle in the 
neighboring patches at rate X/N . When a site x is occupied by i > 1 particles, then it allows birth 
attempts only from particles in the same patch, but an attempt succeeds only with probability 
c(i) (the function c : N — ► [0, 1] represents a self-regulating mechanism as in [5]). We suppose that 
c(0) = 1 (when a site is vacant the first reproduction trial is almost surely successful) and that 
c(j) > c(j + 1) for all j > (the more individuals are present, the less likely it is to produce new 
offspring). Clearly, if c(j) = for all j > iq and c(io) > 0, then in each site there might be at most 
io + 1 particles. 

An equivalent description of the breeding mechanism is that each particle breeds at rate 4> inside 
its patch and at rate A towards each neighboring patch. Newborns choose the target sites inside 
the destination patch at random. The reproduction from a different patch is successful only if the 
target site is vacant, while inside the patch the reproduction is regulated by the function c. 
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Formally, the transition rates at (x, r) for the process rjt are 







1 at rate 



X_ 
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i + 1 at rate 



N 



c(r)(x,r)) Yj r l( x > r i)i f° r 1 ^ 



(2.1) 
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i — 1 at rate 
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In the sequel we refer to the process with these rates as the Patchy Restrained Process or briefly 
PRP. When we need to stress the dependence on the parameters, we write PRP(A, <f>, c, N). We 
observe that many explicit models are recovered as particular cases of this patchy habitat model. 

In particular, by choosing <j> = we obtain a process which does not depend on c and allows 
at most one particle per site. We call this process the contact process on £^, briefly CP. It is 
not difficult to prove that this process has a critical parameter \cp{£%) G [A, Acp(Z d )] (namely 
the population dies out for A below it and survives with positive probability for A above). Indeed 
Acp(^at) < Acp(£i ) = Acp(Z d ) since Xcp{£%) is nonincreasing with respect to N and it is possible 
to prove that Xcp{£%) [ l/2d as N — > oo (see Proposition 13. 7( 1)). 

On the other hand, if we choose c = 5q, we obtain a process that we call biparametric contact 
process (or BCP(A, <j>)) on £fj which coincides with the contact process when = (this is in 
general a process which allows at most one particle per site and has two reproduction rates - the 
inter-patch A and the intra-patch eft). 

The individual recovery process (briefly IRP) introduced in [3] can be seen as a particular case 
of the PRP, namely by setting N = 1 and c = lro,...,/s— li- 
lt is natural to study sequences of PRPs (see Subsection 12.41 and Section [3|) by considering a 
sequence {A^}j G ^ and a corresponding sequence of controlling functions {cj}j g N (we keep A and cj) 
fixed). We are interested in the behavior of the process as i — > oo. We assume that q < Cj + i (the 
regulation gets weaker as i grows). Note that q(j') may be positive for all j, thus there could be 
no a priori bounds on the site carrying capacity. We denote by Coo(i) := lim^^oo c n (i) for all i. 

We conclude this section by listing some particular cases of the PRP which are modifications of 
the IRP. 

2.1. Logistic IRP. If we choose N = 1 and c(i) = max{0, 1 — i/re} (for some kEN \ {0}) then 
we get a cluster of size k at each site in Z d and an interaction between sites which is allowed 
only when the target is empty. This kind of interaction was introduced by Schinazi |llj with a 
cluster recovery clearing mechanism (or mass extinction), and modified by Belhadji and Lanchier 
with individual recoveries in [3] . Note nevertheless that this model differs from the IRP in the fact 
that the breeding inside the cluster becomes increasingly difficult as we approach the full carrying 
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capacity k. For this process r/t, that we call the logistic IRP, the transition rates at site x are: 

— ► 1 at rate A y(z); 

z^x 

( i\ (2.2) 
i — > i + 1 at rate id) I 1 ) , for 1 < i < k — 1; 

V K / 

i — > i — 1 at rate i, for i > 1. 

Clearly the presence of the logistic factor 1 — i/n in the breeding rate may be interpreted as a 
self-regulating mechanism of the population which slows down the reproductions when the patch 
is almost completely exploited. 



2.2. Self-regulating IRP. If N = 1 in general we get a process rjt that we call self-regulating 
IRP. The process has the following rates at site x: 

— > 1 at rate A rj(z); 

z^x 

i — > z + 1 at rate icj)c{i); 

(2.3) 

i — > z — 1 at rate i. 

2.3. Logistic IRP with persistent inter-patches reproduction. Given a general PRP, we may 
identify all the sites in the same patch, meaning that we consider the process £t(x) = SrgKjy Vt ( x ' r )- 
We refer to & as the projection on Z rf of Jft. The projection is a Markov process only when c = (5o- 
In this case it has the following transition rates at site x: 

i^i + 1 at rate Ud) + A ^ ^(z) j ^1 - , for < i < N - 1; 
i — > i — 1 at rate i, for « > 1. 

We note that this process is similar to the IRP of [3] with cluster size equal to N. Nevertheless it has 
two main differences: the interaction between clusters is always active (while in the IRP it vanishes 
once the cluster is nonempty) and it gets more difficult to increase the number of individuals in 
the cluster if the cluster is crowded. 



2.4. The mean field equations for the self-regulating IRP. The usual approach, before 
actually studying the spatial stochastic model, is to derive a non-spatial deterministic version 
called the mean field. We consider the differential evolution equations for the concentrations {tti}i>o 
where U{ represents the fraction of sites having % individuals (clearly ^ m = 1). We compute the 
stationary solutions of the system of equations: a stationary solution with uq < 1 corresponds to 
survival, while no > 1 suggests almost sure extinction. To avoid unnecessary complications, in this 
preliminary study, we take N = 1, that is, there is no difference between patches and sites. 
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Let us discuss briefly the mean field equation for the logistic IRP. The mean field equations are 







1 < i < « - 1, 



u[ = 2n 2 + XuqYa=i iu i -ui(l + (f>(l - £)), 
u' i = {i + l)u l+1 + (i - 1) (1 - <M-i - i (1 + (1 - i) <j>) i 

u' K = -ku k + (k - 1) (l - ^i) <Mc-l- 
Put v! i = 0, for i = 1, . . . , k and sum equations from i to k. Then the solutions, for i = 1, 

i-l 

- — I - I "I • 
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We plug it into the first equation and require u' = to get 



UQ 



=0 (re-i-1)! u 



1 / re— 1 i / . \ 

\ i=l i=l v 



A 

From this we get that uq — ► as k — > oo (thus indicating the possibility of an endemic state for k 
large) when > 1 (<^> > 1 implies exponential convergence, when <p = 1 then no ~ C/X^/k,). 

More precisely, Monotone Convergence Theorem implies that no j min(0, ) as k — > oo. Hence, 
if <f> + A < 1 then no > 1 for all k, while if <j) + A > 1 then no < 1 for k sufficiently large. 

If we consider the self-regulating IRP the mean field equations are 



II n 



Ul 



u[ = 2n 2 + An YluLi iu i ~ui(l+ 4>Cn(l)), 



u; 



(i + l)u i+ i + (i - l)c n (i - l)(j)Ui-i - iui(l + (pc n (i)), 1 < i. 



We can write the solution of the system u\ = 0, for i > 1 as 

ft' 1 

m = — — c n l(i - l)ni, 
i 

where c n \(i) := 131=0 ^(O- We plug it into the first equation and require n' = to get 



1 



n 



Since c n f Cqo as n — > oo then, by Monotone Convergence Theorem, 



1 



> min 



1-0 
A 



and the equality holds if and only if Cqo = 1. In this case, as before for the logistic IRP, <f> + A < 1 
implies no > 1 for all c n , while <j) + A > 1 implies no < 1 for n sufficiently large. 

Thus the mean field model for the self-regulating IRP suggests that, in the spatial case, if 
+ 2d\ < 1 there is extinction for all controlling functions c n (see Theorem 13. If) and, if Cqo = 1, 
survival is implied by <f> + 2dX > 1 provided that n is sufficiently large (we prove a slightly different 
result, see Theorem 13. 5p . 



3. Main results and discussion 



In this section we discuss the effect of the parameters A, (ft, c and N on the behavior of the PRP. 
Part of the arguments are done using monotonicity and coupling with known processes. Indeed 
for any fixed N, the PRP is attractive with respect to A, (ft, c. Moreover we may couple a PRP on 
£fj with another one on £f = Z, d in the following natural way. Consider two PRPs: {rjt}t>o on £% 
(with parameters A, (ft and function c) and {£,t}t>o on ^ d (with parameters X\,(fti and function c\). 
If A > Ai, (ft > 4>ii c > ci then the projection of rj t on Z d (namely Yl r ^K N ^('5 r )) dominates £*(•)■ 

This kind of coupling is the key of the proof of the following result, which states that if the 
breeding parameters A and (ft are sufficiently small then we have almost sure extinction. 

Theorem 3.1. For all functions c, if (ft + 2dX < 1 then there is a.s. extinction. 

By monotonicity it is clear that, if N, c and (ft are fixed then there exists a critical A cr = 
X cr ((ft,c,N) such that given A > A cr the PRP survives with positive probability, while A < A cr 
implies almost sure extinction. By stochastic domination X cr < \cp(£fj) thus if A > \cp{£fj) then 
then PRP survives with positive probability for all choices of (ft and c. 

The following theorem describes the dependence of A cr on (ft and c. Since c{i) can be equal to 
for some i, we identify 1/0 with +oo. 

Theorem 3.2. (1) Given any (ft, c such that Y^'n=o( < l }n Y\i=o c (*)) _1 < +°° then for every A > 

the PRP survives with positive probability, hence X cr ((f),c,N) = 0. 
(2) Given any <ft, c such that Yl^=o ^ Y\a=q c W < +°° then X cr ((f>,c, N) G (0, \cp{£fj))- 

In particular from the previous theorem we deduce that the value of linij_ ) . 00 c(i) (which exists 
by the assumption of monotonicity that we made) tells us whether A cr = or not in almost every 
case. 

Corollary 3.3. Let c(oo) = limj^ 00 c(i). If <p • c(oo) > 1 then X cr ((p,c,N) = 0. If <j> ■ c(oo) < 1 
then X cr ((j),c,N) > 0. As a particular case, for all functions c, 4> < 1 implies X cr > 0. 

In the case <j> ■ c(oo) = 1 the behavior depends on the speed of convergence of c(n) to c(oo). For 
instance if c(n) := 2(^+2)^ anc ^ = 2 then Theorem I3.2f 1) applies and A cr = 0. Unfortunately, due 
to the monotonicity of c, in the case <j> ■ c(oo) = 1, Theorem 13.2( 2) is useless (since the series is 
always divergent). We do not know of any example when A cr (^, c, N) > 0. 

Again by monotonicity, if N, c and A are fixed then there exists a critical <j) cr = (f> C r(X, c, N) such 

that given <p > <f> cr the PRP survives with positive probability, while (ft < <ft cr implies almost sure 

extinction. In general it could happen that (ft cr = oo (almost sure extinction for all (ft - for instance 

when N = 1, c = S and A < A CP (Z d )) or (ftcr — (survival for all positive (ft - for instance when 

A > Xcp{£fj))- The following theorem gives sufficient conditions for cft CT G (0, oo). 
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Theorem 3.4. (1) If A > \ C p{£%) then cp cr (X,c,N) = for all c, N and the PRP survives with 

positive probability when = 0. 
(2) For any A G (0, l/2d) and c such that c(l) > then 4>cr{\ c, N) G (0, +oo). 

Actually to prove <p cr < oo we only need the hypothesis c(l) > while A € (0, l/2d) is needed 
to prove (p cr > 0. One guesses that 4> cr > could be proved under the milder assumption A G 
(0, Acp(£jy)); in order to do so one could adapt the proofs of [3l Theorem 3] and of [TH Theorem 
1(c)]. The main difficulty is to prove the analog of [TJ Theorem 1.7] for the CP on Since 
\qp{£%) I l/2d, extending the interval for A to (0, \cp{£fj)) seems to be a minor improvement at 
least for large values of N. 

Now we want to study the dependence of the behavior of the PRP on the underlying space: note 
that "space" here is described both by N (the horizontal space) and c (the vertical space). Hence 
it is natural to define a sequence of PRPs by means of nondecreasing sequences of functions {cj}j G pj 
and of patch dimensions {AjjjgN- More precisely the ith process is N ^-valued and its transition 
rates are defined as in equation (|2.1[) with c = c, (being d, A, <f> fixed for all i). 

Theorem 3.5. Consider a sequence of PRPs on Efj. with parameters 4> and A and control functions 
Ci. If A > and one of the following conditions holds 

(1) 4>+ 2dX > 1 and N { -> oo; 

(2) 4> > 1 an d Coo = l; 

then there exists i cr = i cr (\,(j),c) such that for all i > i cr there is survival with positive probability. 

Roughly speaking, Theorem 13 . 51 states that if the space is sufficiently large then the PRP survives 
with positive probability provided that the breeding parameters are not too small. In particular 
Theorem 13.51 (1) is a partial converse of Theorem 13.11 In the proof we need to mimic a technique 
introduced in [6] in order to show the "convergence" of the projections on Z d of the PRPs to a 
branching random walk. 

Remark 3.6. Theorems \3.1\ and \3.5\ imply, roughly speaking, that when Ni — > oo, the projection 
of the i-th PRP on 7L d behaves, in the limit i — > oo, as the branching random walk on Z d with 
intra-site breeding rate 4> and inter-site rate X, in the sense that it survives eventually if and only 
if(f> + 2d\ > 1. 

The following proposition deals with the asymptotic behavior of the critical parameters of various 
processes on as Ni — > oo. The results follow from Theorems 13.11 and 13.51 

Proposition 3.7. If N j oo the following hold. 
(1) X CP (£ d Ni ) i l/2d. 
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(2) Consider the biparametric contact process, for all fixed (ft > we have that 

\ BC p{£ d Nl ,4>) imax(i^,o) , 

while for all fixed A > we have that 

<Pbcp (£#.,A) jmax(l-2dA,0). 

In the particular case <j) = A (the intra-patch reproduction has the same rate as the inter-patch 
one), the critical parameter converges to l/(2d+ 1). 

(3) For a sequence of PRPs, 

lim X cr (<j), Ci,Ni) = if 4> > 1, 
lim (/> cr (A,Ci,iVi) = if\> l/2d. 

i— >oo 

We note that Proposition I3.7f 3) gives the limit of X cr as the size of the patch goes to infinity. 
The following proposition deals with the limit of X cr in the case Cqo = 1 (being N and <f> fixed). 
It is a consequence of Corollary 13.31 and Theorem I3.5f 2) and establishes the continuity of A cr with 
respect to c. 

Proposition 3.8. Given N £ N, <fi > 1 and a sequence of functions {cj}j>o such that Cqo = 1, then 

lim X cr (4>,a,N) = X cr (cp, Coo, iV). 

i— »oo 

We recall that the IRP with cluster size k is a particular case of the PRP (taking N = 1 and 
c = l{o,...,re-i})- Thus if X K (4>) is the critical parameter of this IRP, from Proposition 13.81 we get 
the following result on the asymptotic behavior of X K (4>) as k goes to infinity. 

Corollary 3.9. If cf) > 1 then lim^oo A K (</>) = 0. 

We note that, if A°°(</>) is the critical parameter of the IRP with cluster size equal to infinity 
(see [2j for some details on this process) , then as a consequence of Corollary 13.31 if <fi > 1 then 
X°°((f)) = 0. Thus the previous corollary is a result of continuity at infinity for X K (4>) with respect 
to the cluster size. 



4. Proofs 

Proof of Theorem I3.il It is enough to note that the total reproduction rate of each particle is 
bounded from above by 4> + 2dX. Indeed the projection X]fcgK„ Vt(x, k) is dominated by a branching 
random walk on Z d with intra-patch infection rate 4> and inter-patch rate A on each edge. In this 
branching random walk each particle breeds at rate <j) + 2d X and dies at rate 1 . □ 



In order to prove Theorem 13.21 we need the following lemma on discrete time random walks. 
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Lemma 4.1. Let {Z n } n >o be a random walk on N N with the following transition probabilities: 

(ii, .. . ,i N ) — > (h, ... ,ij + 1, ... ,ijv) probability ( f^]y 

(ii,...,ijv) — ► («!,..., - l,...,ijv) with probability jj^jj (4.4) 
(il,...,ijv) — > («'i,...,ijv) with probability (1+ ^ )Af EjliC 1 ~ c fe))- 

fij if SraLo(^ n Eir=o c (*))~ 1 < +°° ^ e random walk is transient. 

(2) The random walk is positive recurrent if and only if X^nLo ^ n llT=o C W < +°°- 

Proof. (1) We note that the projection of the ^-dimensional process on the jth coordinate is a 
birth and death process with forward probabilities <pc(i) / (N (1 + cp)) and backward probabilities 
l/(iV(l+0)) which by [131 Theorem 5.9] is transient if and only if X]^Lo(^ n 117=0 C W) 1 < +oo. 
But transience of at least one projection implies transience of the TV-dimensional process. 

(2) Positive recurrence is equivalent to the existence of a finite invariant measure. It is not difficult 
to prove that 

v(h,...,i N )= n (^n c (*))i f(o,...,o)=i, 

.f->, -0 \ i=0 / 

is a reversible measure. The claim follows noting that 
u(N N ) = Yl v(h,--.,iN) = 



,IN 




□ 



Proof of Theorem \3.2[ (1) Take the PRP with A = 0. Since it cannot leave the original patch, this 
is a (continuous time) random walk on and it is stochastically dominated by the original 
PRP; in this case survival means "avoiding the origin". The survival of the (continuous time) 
random walk is equivalent to the survival of its discrete counterpart (see equation (|4.4p ) which 
is guaranteed by Lemma 14.1( 1) noting that transience implies survival. 

(2) By (1), X cr < Xcp{£ff). We now prove that A cr > 0. Let us start by considering the case 
N = 1 (remember that £f can be identified with Z rf ). For A small, we want to find a subcritical 
branching process which dominates the total number of successful inter-patch reproductions. 
Indeed in this case the PRP cannot survive without leaving the first patch, hence survival is 
equivalent to the positive probability of having an infinite number of successful inter-patch 
infections. Once a particle appears at x, the total number of particles living at x up to the next 
time when the colony at x has no individuals is, according to Lemma 14.11 a positive recurrent 
random walk (let us call tq the number of steps before reaching 0) with the following rates: 

i — > i + 1 with rate i(j)c(i) 

i — > i — 1 with rate 1 

i — > i with rate i(j){l — c(i)). 
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Note that the exponential clock of the transition (including the loop) has parameter i(l + </>). 

The number of inter-patch reproduction trials made from the site x depends on the number 
of particles at x: indeed, when the number of particles at x is i, the total reproduction rate 
towards the neighboring patches is iX. Hence the number of inter-patch reproduction trials 
between two transitions does not depend on i, since it is the number of times a Poisson clock 
with rate iX rings before the transition clock with rate i(l + (j>) does. It is clear that the 
total number of reproduction trials (before the progeny of the original particle at x dies out) 
from x to the neighbors is YlkLi S?=i ^fc(i) where {Yfc}fc>i is an i.i.d. sequence of random 
2<i-dimensional vectors; to be precise, if m is the number of reproduction trials on the neighbor 
% then 

(E-=i^)!(A/(l + 0))^ 



F(Y k = (m,...,n 2d )) = , ^, 



(l + 2dX/{l + 4>)) l +^n l 



n 



(see [7] for details). Thus 



as A —* 0. By choosing A < (1 + <p) / (2(IK(tq)) the total number of successful trials (which 
cannot exceed the total number of trials) is dominated by a (subcritical) branching process 
with expected number of offsprings 2dj^E(ro) < 1- 

Let us take N > 1: we cannot repeat exactly the same argument since generations inside a 
patch are not independent, due to the presence of the controlling function c. Hence we couple 
with a suitable process. The coupled process has the following rules: 

a. every inter-patch reproduction is successful (i.e. we do not use c); 

b. we define the particle born from an inter-patch reproduction as the ancestor of its descen- 
dance inside its patch; 

c. inside a patch particles compete for resources if and only if they have the same ancestor 
(i.e. the function c applies inside the same progeny). 

Clearly this process dominates the original PRP and progenies of different ancestors are inde- 
pendent. 

The number of generations of the progeny of a particle in the dominating process is equal to 
the absorption time tq of a randow walk with rates 

(•&!, . . . , ijy) — ► (ii, ■ ■ ■ , ij + I, • • • , ijv) with rate c/>c(ij) 
(ii, . . . , ijv) — ► {ix-, ... ,ij — 1, . . . , in) with rate 1 
(h, ...,in) — ► (h, ■ ■ ■ , *jv) with rate (pJ^jLii 1 ~ c (*i))- 
As before, the behavior of this random walk is equivalent to the behavior of its discrete-time 
counterpart which, according to Lemma |4.H is positive recurrent (thus E(to) < +oo). Hence 

survival is equivalent to the positive probability of having an infinite number of progenies. But 

ii 



the expected number of successful inter-patch reproductions is equal to 2<iY^E(ro). Choosing 
A small enough yields the conclusion. 

□ 

Proof of Theorem\3.4\ (1) It follows from the fact that for all > the PRP stochastically domi- 



nates the CP on £fj. 

(2) The hypothesis c(l) > yields (j) cr < oo. Indeed let </> C r(A, 2) > be the critical threshold in [31 
Theorem 3] and choose (ft > <^cr(^> 2)/c(l). Let rjt : 7L d — ► {0, 1, 2} be the IRP with parameters 
A = A, (ft = 0c(l) and k = 2. By these choices of the parameters, rjt is supercritical. On the 
other hand, the projection = SigK n Vt(x,i) dominates rjt and this proves that the PRP 
survives when (ft > (ft cr (X,2)/c(l). 

To show that (ft cr > in the case where A < l/2d, note that by Theorem 13.11 if we choose a 
positive (ft smaller than 1 — 2cfA then the population dies out. 

□ 

Before proving Theorem 13.51 we prove an equivalence between conditions on the sequence of 
parameters. 



Lemma 4.2. The condition 



is equivalent to 



mm 



Ni -> oo V Coo = 1 (4.5) 



E °i(fW > 1 (is i > oo, VMS N. (4.6) 



/€N K iV : £,/(j)=M l\i ^ 

Proof. We first prove that (|4,5p implies (|4.6p . Since c, = 5q, if Ni ^ oo then 

1 Ni 

min — 5^c i (/(i)) = l-M/iV ij 

thus condition (|4.6p holds. For general Cj the same follows from a > 5q. If Ni is bounded then 
Ni = N eventually for some JV £ N (since the sequence is nondecreasing) . Hence for any / such 
that J2j fti) = M we have ^ Y^Li c *(/(i)) ~* SjLi Coo(/(j)) which is equal to 1 if = 1. 

Suppose now that (j4.5j) does not hold, that is JVj is bounded and Coo(M) < 1 for some M G N. 
Choose / such that /(l) = M and elsewhere. Thus, eventually 

K min -L £ a(m) < = 1 ~ 1 " C " (M) < 1, 

/eN K iv : £, f(j)=M Nij^ N ^ AT 

which proves that (14. 6p does not hold. □ 

Proof of Theorem 13.51 (1) The PRPs of the sequence dominate the ones where the control functions 
are equal to 5q for all i. Hence it is enough to prove the statement for the latter case. 
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The projection of this process on 'L d , namely Ct( x ) := YlreK r lt( x i r )-> nas transition rates 



at x E Z d 



+ 1 at rate fj</> + A £ ${y) j (l " ^) 



j — ► j — 1 at rate j. 



NiJ ' (4.7) 
This process dominates eventually (as Ni — > oo) the process {£,t}t>o with rates 

(4.8) 



j -► j + 1 at rate I j'0 + A ^ J (1 - e)l [0 ,jv] (i); 



j i - 1 at rate j; 

for all iV > and e > (it suffices that N { > N/e). Choose e > such that (<£+2cZA)(l-e) > 1. 
We will prove now, following [6j, that for all sufficiently large N, the process {£,t}t>o survives 
with positive probability. The strategy is to study the branching random walk {£,t}t>o with 
rates 



j j + 1 at rate lj(j> + A ^ J 



(4.9) 



j — ► J — 1 at rate j, 



starting with one particle at the origin G Z d , where 4> = (1 - e)0 and A = (1 - e)A. Note that 
the process {^t}t>o can be seen as a truncation (at height N) of {Ct}t>o- 

It is not difficult to verify that E <5 °(^) satisfies the following differential equation (see [6, 
Section 5], where we did the same with the classical BRW): 

j t ^ 5o (Ux)) = -E s °(t t (x)) + d>E s °(Z t (x)) + \J2® S °(Uy)), (4-10) 
whose solution is 

Ik \n—k+n 

^ k \0,x)- 

n=0 fc=0 

where ^ n ' k \§, x) is the number of paths from to x of length n and k loops. Moreover, if 
\x\ Z d = 1 then, taking t = n large enough, 

E 5 ° (&(*))> T (0, x) e -t = y (n,k )(Q } J> k \ n ~ k n n $ + 2d\) n 

vs> v " - v ; n! ^ r v J ((j) + 2dX) n n! 

fc — k — 

™ (^ + 2 ^) n yyn.fc) f0 x) ^ n ~ k > (J + MXT d/2 
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,../r \n—k-in 

E 5o (6(*)) = £ E ^ fc) (0,x)^A_± e -*, (4.11) 



for some C = C(A, 6) > 0, since Y2=l (°> x)6 k X n ' k /(6 + 2d\) n is the probability of being 
in x (\x\%d = 1) after n steps for a discrete-time random walk with transition probabilities 



p(x,y) 



A f 

if x ~ y, 



+ 2dX 

6 

ifx = y, 



b + 2dX 

(see [121 Corollary 13.11] for the asymptotic estimates of the n-step probabilities). Hence there 
exists t = n sufficiently large such that E °(^(a;)) > 1 (this result is analogous to Lemmas 5.2 
and 5.3 of [6j). Using the same arguments as in Lemma 5.4, Remark 5.5 and Theorem 5.6 of 
[6J, it is possible to prove that the process {Ct}t>o survives when N is sufficiently large and this 
implies the survival of {£t}t>o when i is sufficiently large. 
(2) Let k c (X, (1 + 6)/2) as in O Theorem 5]. By Lemma 14.21 we may pick i cr such that for all 

1 Ni 1 + 6 

min Y,<*(fU))>-£r, VM<«c + l. 

In this case, if in a patch there are at most k c + 1 particles then the effective reproduction rate 
for each particle within the patch is at least (1 + 6)/2, which implies that the PRP dominates 
the IRP with k = k c + 1 and reproduction rates A and (1 + 6)/2 respectively. 

□ 



5. Open questions 

In the phase diagram the answer to some questions remains unknown. For instance, we do not 
know if X cr (6,c, N) = when 6 ■ c(oo) = 1 (we know there are examples where this happens). 
Moreover we do not know if 6 cr (X,c,N) > when A G (l/2d, Xcp(Sff))- One interesting question 
is whether a nontrivial invariant measure exists when there is survival with positive probability. 
The answer is positive if c(i) = eventually (in this case the configuration space is compact and the 
invariant measure is obtained as the limit of the process starting from the maximal configuration) . 
In the noncompact case one could try to mimic the technique employed in [5]. 

On the other hand, one may consider other models which, although similar to those considered 
here, cannot be obtained as particular cases of the PRP. One of particular interest is a self-regulating 
IRP where the inter-cluster infection is always active. In each site of Z d there is a (possibly infinite) 
cluster and, given a regulating function c, the transition rates at site x are: 

i — > i + 1 at rate i^i6 + A f/(z)^ c(i) , for < i; 
i — > i — 1 at rate i, for % > 1. 
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In the present work, we dealt with the case c(i) = max(0, 1 — i/N) (see Subsection 12. 3p . Further 
investigations may be made on the processe where we add a rate 5 for the transition i — ► - the 
so-called catastrophe or mass extinction. 
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